Scenarios of heterogeneous nucleation and growth studied by cell dynamics simulation 
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The dynamics of phase transformation due to homogeneous nucleation has long been analyzed us- 
ing the classic Kolmogorov- Johnson- Mehl-Avrami (KJMA) theory. However, the dynamics of phase 
transformation due to heterogeneous nucleation has not been studied systematically even though 
it is vitally important technologically. In this report, we study the dynamics of heterogeneous 
nucleation theoretically and systematically using the phenomenological time-dependent Ginzburg- 
Landau (TDGL)-type model combined with the cell dynamics method. In this study we focus on the 
dynamics of phase transformation when the material is sandwiched by two supporting substrates. 
This model is supposed to simulates phase change storage media. Since both homogeneous and het- 
erogeneous nucleation can occur simultaneously, we predict a few scenarios of phase transformation 
including: homogeneous-nucleation regime, heterogeneous-nucleation regime, and the homogeneous- 
heterogeneous coexistence regime. These predictions are directly confirmed by numerical simulation 
using the TDGL model. The outcome of the study was that the KJMA formula has limited use 
when heterogeneous nucleation exist, but it could still give some information about the microscopic 
mechanism of phase transformation at various stages during phase transformation. 

PACS numbers: 64.60.Qb, 68.18.Jk, 81.10.Aj 
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I. INTRODUCTION 



The dynamics of phase transformation by nucleation 
and subsequent growth of a nucleus of a stable phase is a 
very old problem, which has been studied for more than 
half a century from a fundamental point of 

view as well as from technological interests. However, it 
plays a vital role today, for example, in information tech- 
nology (IT) of phase-change optical data storage media. 
For example, rewritable (RW) compact disks and digital 
video disks (DVD) consist of Te based alloys. The data 
unit "bit" is switched between and 1 by laser heating 
between the amorphous and the crystalline state i^, [S] . 
The crystallization occurs through the heterogeneous nu- 
cleation at the surface of the disk. Therefore, the data 
transfer rate of such a phase change storage media de- 
pends critically on the time scale of amorphous to crys- 
tallinephase transformation by the heterogeneous nucle- 
ation d, Q. 

The dynamics of phase transformation is usually 
described by the Kolmogorov- Johnson-Mehl-Avrami 
(KJMA) kinetics rule 1, 2, 3, ,4,^]. According to this rule 
the fraction of transformed volume follows a straight line 
in the so-called Avrami-plot and the dynamics is char- 
acterized by the slope called the Avrami exponent. The 
time scale of phase transformation is determined from 
the incubation time that is the waiting time before the 
start of the phase transition, and the growth time that 
is the time necessary to accomplish the phase transfor- 
mation. The incubation time depends on the transient 
nucleation rate 0, while the growth time depends both 



'Electronic address: ' i wamat su@ph . ns . musashi-tech. ac . j p | 



on the steady-state nucleation rate and the growth rate 
of a nucleus. 

Theoretically, the validity of the KJMA kinetics is well 
understood [l|, |3, S 3 and tested [13, UH for homoge- 
neous nucleation, but it is not always obvious for hetero- 
geneous nucleation [l| . Various modeling and analytical 
calculations seem to have limited use [IJl- Therefore, 
theoretical modeling and a computer simulation method 
that can handle both homogeneous and heterogeneous 
nucleation on the same footing are highly desirable. 

In the previous paper [Til [ij] , we adopted the cell dy- 
namics method [lj| to study KJMA kinetics of phase 
transformation within the framework of time-dependent 
Ginzburg-Landau [3] or phase field [lB| model. The 
cell dynamics method can successfully describe the clas- 
sical picture of homogeneous nucleation as well as the 
interface-limited growth [l3| of nuclei. Therefore, the 
whole process of KJMA kinetics [H of homogeneous nu- 
cleation and growth can be simulated on the same footing 
using cell dynamics. Then, the time scale of transforma- 
tion that consists of growth time and incubation time can 
be studied. The same story should be effective for hetero- 
geneous nucleation and growth as well. Similar studies 
of heterogeneous nucleation usin g th e phase field model 
have appeared very recently [isl Il9l. [20 | , but they paid 
more attention to the early stage of grwoth by spreading. 

In this paper, we will use the cell dynamics method 
with thermal noise to study heterogeneous nucleation 
and growth under various conditions and settings in the 
TDGL model in a unified manner. In Section II we 
present a short review of the classical picture of nucle- 
ation and growth. We also discuss a few expected sce- 
narios of heterogeneous nucleation present on the surface. 
In Section III, we will present the cell dynamics method 
with thermal noise for studying heterogeneous nucleation 
and growth. The surface field induces heterogeneous nu- 
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cleation. Section IV is devoted to the results of numeri- 
cal simulations. There, we pay attention to the Avrami 
exponent, the incubation time, and the nucleation rate 
which play dominant roles in the transformation time 
scale. Finally Section V is devoted to the conclusion. 



II. DYNAMICS OF HETEROGENEOUS 
NUCLEATION AND GROWTH 

The time evolution of the volume fraction /homo of 
transformed volume for homogeneous nucleation, for ex- 
ample in two-dimensional infinite system, predicted h 
KJMA (Kolmogorov-Johnson-Mehl- Avrami) [U, 0, S " 
0j theory is as follows: 



/ = 1 - exp ( -^/s^s {t - Uncf 



which can be generalized as 



/ho 



1 — exp 



(1) 



(2) 



where m is the so-called Avrami exponent, ti^c the incu- 
bation time, and Tgr is the growth time, which is deter- 
mined from the steady state nucleation rate Is and the 
steady state interfacial velocity Vs through 
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(3) 



Then, the double logarithm of ([2]) gives a straight line 

logio - /homo)) = m logio {t - tine) + const. (4) 

between log (— ln(l — /homo)) and log {t — t-mc) which is 
the so-called Avrami plot. The slope gives the Avrami 
exponent m. 

According to the theory of Shncidman and Wein- 
berg , if we take into account the time-dependent tran- 
sient nucleation rate for Is as well as the size-dependent 
interfacial velocity for Vs and the finite size of the critical 
radius iZ* of nucleus, the incubation time tmc is given by: 



kBT 



(5) 



where t ~ R*/vs is the time lag due to the finite size i?, 
of the critical nucleus, is the energy barrier to form 
the critical nucleus and T is the absolute temperature 

In classical nucleation theory P, [U, HI, HI] i *he steady- 
state homogeneous nucleation rate Is , that is the number 
of critical nuclei which appear per unit time and unit 
volume, is usually given by the activation form 



/.occxp(-i^). 



(6) 



where W, is the nucleation barrier of critical nucleus that 
appeared in 

In order to derive KJMA formula ^ we have used the 
following assumption |9(|: 



(a) the system is infinite and no boundary effect exists. 

(b) the nucleation is homogeneous. Therefore, the nu- 
clei form at random sites. 

(c) the steady state interfacial velocity Vg is constant 
or weakly (power-law) time-dependence. 

(d) the steady state nucleation rate Is is constant or 
weakly (power-law) time-dependent. 

In an ideal case when the steady state velocity Vg and 
the steady state nucleation rate are strictly constant, the 
Avrami exponent to in ^ is given by 



m — d + 1 



(7) 



for c?-dimensional system. Then to = 3 for d = 2 di- 
mensional system. Incidentally, when the fixed num- 
ber of nuclei exist initially and continuous nucleation is 
suppressed (no nuclei are added), which is called site- 
saturation case we have 



= d 



(8) 



Therefore, the Avrami exponent becomes smaller in the 
site-saturation case. 

The first and the second conditions (a) and (b) are 
violated, however, for heterogeneous nucleation. Fur- 
thermore, the energy barrier is lowered by the 
amount 



with 



fie) < 1 



(9) 



(10) 



for the circular meniscus, where 9 is the contact an- 
gle [2§|. For the complete or partial wetting condi- 
tion 9 = 0, the wetting layer spread over the surface 
which means that heterogeneous nucleation occurs with- 
out crossing the energy barrier since W^,f{9) = from 
pop . Furthermore the incubation time becomes negative 
from (pS)) . Even in the partial wetting [2^, [2^ condition 
< 9 < tt/2, the energy barrier is lowered W^f{9) < W* 
and heterogeneous nucleation can occur more easily than 
homogeneous nucleation. The incubation time for het- 
erogeneous nucleation also becomes shorter than for the 
homogeneous nucleation [24| . 

When heterogeneous nucleation can occur near the 
substrate, we may expect coexistence of heterogeneous 
nucleation near the substrate and of homogeneous nucle- 
ation in the bulk. Then the time evolution of the volume 
fraction / of transformed volume would be given approx- 
imately by the superposition of the two contributions: 



/ — /ho 



/i 



hot cro 



(11) 



where /hotcro is the transformed volume by homogeneous 
nucleation ^ and /hctcro is the one by heterogeneous 
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nucleation. The transformed volume fraction increases 
due to the simultaneous occurrence of both homogeneous 
and heterogeneous nucleation. 

The latter contribution can be estimated by consider- 
ing when the homogeneous nucleation is completely sup- 
pressed. Then, we may expect that nucleation occurs 
only at the substrate as heterogeneous nucleation and 
new phase growth from substrate as a traveling wave 
with front velocity Vs- Then, the transformed volume 
is simply proportional to the time: 



yhctcio Vgt 



(12) 



since the front velocity Vg depends only on the under 
cooling e of the material from pO|) and does not depend 
on the existence of the substrate. In this case, the usual 
KJMA formula ^ cannot be used. 

Summarizing, we may expect a few phase transforma- 
tion scenarios: 

(a) The homogeneous nucleation regime where hetero- 
geneous nucleation can be neglected. The phase 
transformation is governed by homogeneous nucle- 
ation in the bulk. This situation occurs when the 
substrate is no- wet {9 = tt). 

(b) The heterogeneous nucleation regime where homo- 
geneous nucleation can be neglected. This can oc- 
cur when the temperature is low or the under- 
cooling e is low and the substrate is wet {6 < tt) 
because the homogeneous nucleation rate will be 
low and would not be visible. 

(c) The coexistence regime where both homogeneous 
and heterogeneous nucleation can exist. This can 
occur when the temperature is high or the under- 
cooling e is large and the substrate is wet {9 < tt), 
then the homogeneous and the heterogeneous nu- 
cleation rate as well as the interfacial velocity could 
be high. 

In those cases, the popular KJMA kinetic rule given by 
^ is not applicable. We have to analyze the transformed 
volume fraction / directly. 



III. CELL DYNAMICS SIMULATION FOR 
HETEROGENEOUS NUCLEATION AND 
GROWTH 

We use the standard time-dependent Ginzburg Lan- 
dau (TDGL) [ill model to study the dynamics of het- 
erogeneous nucleation. The cell dynamics simulation 
method is used to solve the TDGL equation numerically. 
This method digitizes the time and space of standard 
TDGL [l^ evolution equation 



'dt 



5T_ 



(13) 



where 5 denotes the functional differentiation, ■0 is the 
non- conserved order parameter, and T is the free energy 
functional. This free energy is written as the square- 
gradient form 



= \j [D{Vi^f -f /i(^)] dr. 



(14) 



The local part h{il)) of the free energy functional deter- 
mines the bulk phase diagram and the value of the order 
parameter in equilibrium phases. 

This TDGL equation (fT3l) is loosely transformed into 
a space-time discrete cell dynamics equation [14]: 



iP{t + l,n)=F[,p{t,n)]+^{t,n), 



(15) 



where the time t is discrete and an integer, and the space 
is also discrete and is expressed by the integral site index 



The thermal noise ^(i) is added to ([13]) in (fT5)) 
which is related to the absolute temperature T from the 
fluctuation-dissipation theorem as 



(16) 



In this paper, we will use a uniform random number rang- 
ing from to -1-^0 Then, the parameter is 
proportional to the absolute temperature: 



(17) 



or the temperature T is included through the thermal 
noise ^o- 

Since the cell dynamics method is invented not to simu- 
late the mathematical TDGL partial differential equation 
([13]) accurately but to simulate and describe the global 
dynamics directly |14| . it certainly cannot simulate the 
high-frequency fluctuation in space and time. Therefore 
the thermal fluctuation eq. (jl7p cannot take into account 
the short-time fluctuation. However, such a fluctuation 
will be unimportant since the thermal fluctuation plays 
role mainly in the birth of nucleus and will play secondary 
role in time-evolution of nucleus. 

The mapping F is given by 

F[^{t, n)] = ~f[^{t, n)) + D [« 0(t, n) > -iP{t, n)] , 

(18) 

where 



(19) 



is the derivative of the local part /i(V') in ([14)) . and the 
deflnition of <C ■ • • ^ for the two-dimensional square grid 
is given by |1J| 

« V(i,") »= ^ ^ ^(i,i) + ^ Yl ^'(i,^, (20) 
i=nn i=nnn 

where "nn" denotes nearest neighbors and "nnn" next- 
nearest neighbors. In this report, we consider only the 



4 



two-dimensional square grid to avoid heavy numerical 
work. 

We use the map function f{ip) directly obtained from 
the free energy /i(V') from (HH) instead of the stan- 

dard tanh form originally used by Oono and Puri [l^ . 
which is essential for studying the subtle nature of nu- 
cleation and growth when one phase is metastable and 
another is stable. 

In order to study heterogeneous nucleation, we con- 
sider the situation where the upper and lower side of 
material is sandwiched by substrates. This situation sim- 
ulate heterogeneous nucleation observed in Te based al- 
loy of rewritable (RW) compact disks and digital video 
disks (DVD) [ESi]. The effect of the substrate, which 
induces heterogeneous nucleation, is represented by the 
surface field, which has been used to study the wet- 
ting [2^ [2^ and heterogeneous nucleation 27]. Here, 
we consider a x — z two-dimensional system, and assume 
that the material is sandwiched by two lines parallel to 
the a;-axis at 2; = and z = d. 

The local part of the free energy /i(V') we use [ill, E^l 
consists of two parts: 



h{il;)^ho{^) + hs{^) 
where the bulk free energy hQ{ip) is given by 

>3 



T 



(21) 



(22) 



This free energy is shown in Fig. [1] where one phase at 
ipm = is metastable while another phase at V's = 1 is 
stable. The free energy difference Ah between the stable 
phase and the metastable phase is determined from the 
parameter e: 

Ah = h{^jj,n = 0) - hi-^, = 1) = ^. (23) 

Therefore, e represents the undercooling which could be 
measured from melting point. We will use the terminol- 
ogy undercooling to represent e. The metastable phase 
at tpm — becomes unstable when e = 1/3 = 0.33 . . . , 
which defines the spinodal. 

The surface free energy hs{ip) is given by [l^H^I 



h.W = - hiHz) + JuS{z - d)) ip 



(24) 



where 7; denotes the short-ranged surface field of the 
lower substrate at 2; = and 7^ denotes one of the up- 
per substrate ai z — d. S{x) is the usual delta function. 
At the surface of the substrate free energy is lowered by 
the amount "fiipi at the lower substrate and ^ui^u at the 
upper substrate, where ipi and ipu are the order parame- 
ters in contact with the upper and the lower substrates. 
Therefore, the substrate prefers to be wetted by the new 
stable phase with ipg = 1 when 7;_u > 0, and the nu- 
cleation of the stable phase occurs preferentially at the 
surface of the substrate. Then, heterogeneous nucleation 
could be promoted at the substrate. 
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FIG. 1: Model double-well free energy /io(^/') defined by ([22|) 
that can realize two-phase coexistence when e = 0. The un- 
dercooling parameter e determines the free energy difference 
Ah between the depth of two wells. The phase with ^prn = is 
metastable while the one with ifjs ~ 1 is stable. The spinodal 
occurs when e = 1/3 = 0.33 . . . 



Due to the presence of the boundary at z = and 
z = d, the TDGL equation (fTS]) should be augmented 
with the boundary condition: 



-'1 



dip 
dz 



-11 



2D 



In 



(25) 



where subscript I means that the derivative should be 
calculated at the surface of the lower substrates at z = 
and u means the upper substrate surface at z = d. 

A similar boundary condition with 7 = was used by 
Castro [31 and later by Granasy [l^] , which corresponds 
to the free boundary condition. It is also possible to fix 
the order parameter, for example, by setting -0; = 1 at the 
boundary. This strategy was also adapted by several re- 
searchers to study the heterogeneous nucleation ^20", '2^ 
and wetting [2§|. We use the boundary condition pS)) 
because it is directly related to the contact angle of the 
popular theory of wetting [2^ [2^. It should be noted, 
however, this contact angle does not mean naive geomet- 
rical angle, rather the angle represents wettability de- 
rived from the thermodynamic consideration. Therefore, 
the boundary condition used by Semoroz [sO ] is merely a 
phenomenological condition and does not represent any 
thermodynamic meaning. 

In cell dynamics code, (|25|) is transformed into a dis- 
crete form. If we use symbolically the index n = as 
the bottom of surface at 2; = and n = N as the top 
of the surface at z = A^, the boundary conditions for the 
phase-field iplt, x, z) are given by 



iP{t,x,0) = ^{t,x,l) + ^, 
ip{t,x,N) = ip{t,x,N - 1) + 



7u 
2L>' 



(26) 
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FIG. 2: A typical evolution pattern of the phase transformation for 512x512 system sandwiched by a neutral substrates with 
7u = at the top and a non-wet substrate with 7; = —0.4 at the bottom calculated from our cell dynamics code for a medium 
undercooling e = 0.2 and at a lower temperature I/^q = 60. Nucleation occurs predominantly within the bulk and the nuclei 
grow isotropically preserving an almost circular shape. At a later stage {t — 200) a small amount of metastable materials are 
left as bubbles of metastable phase (black) attached to the lower substrate. 



which takes into account the effect of the short-range field 
7„.; of the surface. In this case, however, the function / 
in the map function F in (jl8[) is given f{ip) = dhQ{ip) / dip 
of the bulk because the effect of the surface field is short- 
ranged. A similar condition was used by Marko [31] in 
his cell dynamics model to study the surface effect on 
the dynamics of the phase transformation of a conserved 
system. 

By changing the magnitude of the interfacial energy 
7„ and 7;, we can change the heterogeneous nucleation 
rate on the substrate. It should be noted, however, that 
the growth velocity of the nucleus, even if they are het- 
erogeneously nucleated at the substrate, depends only on 
the undercooling e of the bulk The substrate acts 

only as the seed of the nuclei, therefore it only affects the 
nucleation rate at the substrate and does not influence 
the growth of the nuclei. 



which can be calculated from the interfacial profile cal- 
culated from the Ginzburg-Landau equation 6!F/S^ — 
for the order parameter ip at the two-phase coexistence 
(Aft, = 0) [231 . Then the e-dependence of the surface 
tension a could be ignored, and the energy barrier W^, 
in (|27p is mainly determined from the magnitude of un- 
dercooling e. The energy barrier tV* is high when the 
undercooling e is low. 




-0.4 



60 



IV. NUMERICAL RESULTS AND DISCUSSION 

We used the cell dynamics code developed previ- 
ously [Tl| . and simulated phase transformation when 
the effects of substrate and of the heterogeneous nucle- 
ation exist. In our TDGL model (HH) with (gl]), the en- 
ergy barrier of the homogeneous nucleation of two- 
dimensional circular nucleus is given by 



Ah 



1 

oc — 
e 



with its critical radius i?* given by 

a 1 



Ah 



(27) 



(28) 



where Ah is the free energy difference (f23| between the 
metastable and the stable phase (Fig. [1]). The interfacial 
energy a between the stable and the metastable phase is 
given by 



(29) 
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FIG. 3: The time evolution of the volume fraction / average 
over 50 samples for 512 x 512 system when the material sand- 
wiched by a neutral substrates with 7^ = at the top and 
a non-wet substrate with 7; — —0.4 at the bottom (broken 
line) compared with the values when the system is infinite 
(homogenous nucleation). Standard deviations are not shown 
because they are too small to be visible in this scale. 

In our TDGL model, it has been shown that the nu- 
cleus grows by the interface-limited growth. The inter- 
facial front velocity Vg is almost constant and is roughly 
given by [ll,|I3,[3| 



(30) 



using the undercooling e. The deeper the undercooling, 
the higher the front velocity becomes. The interfacial 
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FIG. 4: The Avrami plot corresponding to Fig. |3l Again two 
curves are almost all in straight lines and the slope (Avrami 
exponent m) are almost the same. The Avrami exponents 
m deduced from the linear portion of this Avrami plot are 
m=4:A6 for 7=-0.4. This value is very close to m = 4.45 for 
homogeneous nucleation. 

velocity Vs does not depend on the temperature. There- 
fore, the temperature dependence of growth time Tgr in 
comes solely from the temperature dependence of nu- 
cleation rate Is given by Furthermore, the growth 
velocity of nuclei, which nucleate in the bulk and sub- 
strate, should be the same. 

The contact angle is connected to the surface energy 
7 through the Yound-Laplace equation [12, 

COS0=2^(V.,-Vrn) = ^ (31) 

a a 

as the surface free energies are given by and 7;,„V'm 

respectively at the lower and the upper substrates when 
the substrates are covered by the stable phase ips = ^ 
and the metastable phase ip^n = 0. 

Thermal noise is introduced as a uniform random num- 
ber r picked up from — < Co- The system size is 
fixed to 512 x 512 since a larger system with 1024 x 1024 
was shown to give almost the same results [ll| . Through- 
out this study, we set the value D = 1/2. Then the sur- 
face tension a between stable and metastable phase is 
cr= 1/24- 0.042. 

A. Homogeneous nucleation regime 

When the substrate is non-wet 7 < 0, heterogeneous 
nucleation is completely suppressed. Then, phase trans- 
formation occurs through homogeneous nucleation and 
growth. 

Figure[2]shows a typical pattern of evolution of a stable 
new phase (white) in the meta-stable old phase (black) 
when the material is sandwiched by the supporting sub- 
strate with a non-wet substrate with 7; = —0.4 on the 
bottom. The top surface is supported by the substrate 
with 7„ = 0, which means a free boundary condition is 



imposed [T^. The contact angle 9 at the bottom sub- 
strate is determined from cos6' = —0.4/0.042 ~ —10, 
which gives the completely drying contact angle 9 = n. 
Heterogeneous nucleation or wetting at the substrate is 
completely suppressed. 

The contact angle at the neutral top substrate is given 
hy 9 = n/2. We consider medium undercooling e = 0.2 
and low temperature 1/T oc l/^g = 60. The periodic 
boundary condition is imposed along the horizontal x- 
axis. Similar to the homogeneous nucleation considered 
previously [ll|, the nucleation occurs inside the bulk. 
The shape of the nucleus is almost circular and will grow 
isotropically with a seemingly constant interfacial veloc- 
ity Vs ■ Therefore, the classic picture of homogenous nu- 
cleation and growth seems valid in this honiogeneous- 
nucleation regime. 

Figure [3] shows the time evolution of the volume frac- 
tion / of transformed volume for two-dimensional 512 x 
512 system with 7; = —0.4 averaged over 50 samples. 
Since nucleation occurs uniformly as shown in Fig. [51 the 
sample-dependence of / is small. The standard devia- 
tions for / at each time step is too small to be visible in 
Fig. O The curve is very close to that of homogeneous 
nucleation in the infinite system, which is also averaged 
over 50 samples. The curve for 7; = —0.4, however, can- 
not be fitted to theoretical curve ([2]) by the least-square 
fitting of three parameters, Tgr, tine and m in principle 
because system is not homogeneous any more. However 
it could be well represented by the KJMA formula ^ 
since the evolution of / does not differ much from that 
of homogeneous nucleation. 

In fact, our simulation data in Fig. [3] can be converted 
to the predicted KJMA linear relation ^ as shown in 
Fig. m In order to convert the volume fraction / to the 
Avrami plot, we have to specify the incubation time ii„c 
in (g]). Since we cannot fit the KJMA formula (g]) di- 
rectly to the simulation data for /, we cannot use least- 
square fitting to deduce tjnc- Instead, we have directly 
deduced incubation time from the averaged data of / 
shown in Fig. [31 We set tine = ^o where to is the time step 
when the volume fraction / start to increase from 0. The 
transformed Avrami plot shows linear relation predicated 
from ^ which means that the classical KJMA formula 
^ seems to be correct in this homogeneous nucleation 
regime. The exponents m deduced from this Avrami plot 
are almost the same for different 7 and are m =4.36, 4.26 
and 3.86 for 7 = 0.0, -0.02 and -0.4 respectively, which 
are very close to the m =4.45 for the homogeneous nu- 
cleation. These values, however, are larger than the ideal 
value m = 3 calculated from ([7]) when d — 2. 



B. Heterogeneous nucleation regime 

In contrast to the homogeneous nucleation regime in 
Fig. [21 a new nuclear embryo preferentially appears near 
the bottom substrates and the new phase starts to grow 
from the bottom substrates in the heterogeneous nucle- 



t=l2(Xl 



t=6lXlQ 



FIG. 5: The same as Fig. [2] but with a weakly wet substrate with 7; = +0.07 at the bottom for the low undercooling e = 0.1 
and at a medium temperature 1/Co = 40. Homogeneous nucleation is suppressed and the nucleation occurs preferentially from 
the lower wet substrate. Only a very small number of nuclei formed at the substrate because of weak surface field, which leads 
to the very strong sample dependence of the evolution of the volume fraction /. 



ation regime. When the surface field 7; is weak, only a 
small number of nuclei is formed at the lower substrate 
and grow as shown in Fig.O Initially we observe semicir- 
cular nucleus nucleated on the bottom substrate. Even- 
tually they grow with constant velocity Vs given by pO|l 
and start to coalesce. Then the evolution of the volume 
fraction / shows large sample-dependence because it de- 
pends strongly on when and how much the nuclei are 
formed. 

Since there is large scattering of data for the evolution 
of volume fraction /, we show in Fig. [6] the time evolution 
of / averaged over 50 samples. There are large standard 
deviations for 7; = 0.05 and 0.07 which could be as large 
as the average / itself. This is because there is large 
scattering in the time when the nucleation starts and in 
the number of nuclei formed. It should be noted that 
the complete wetting condition is satisfied if 7; > 0.042 
as the surface tension of the nucleus is given by cr ~ 
0.042. As we increase the strength 7; of the surface field, 
the scattering of data decreases as the nucleation occurs 
almost uniformly on the substrate as shown in Fig. [T] 
The standard deviation for 7/ = 0.1 and 0.4 becomes less 
than 10% of / itself and have not been shown in Fig. [71 

When the surface field is strong, heterogeneous nucle- 
ation occurs easily, and there are numerous nuclei formed 
at the substrate. Initially we observe a semicircular nu- 
cleus nucleated on the bottom substrate. They soon coa- 
lesce and will grow with constant velocity. Finally the lin- 
ear horizontal moving front of transformed phase appears 
and moves upward with constant velocity Vs (Fig. [7]). 

Plane- wave like moving front reflects in the non-KJMA 
behavior of the volume fraction / of transformed vol- 
ume for 7/ = 0.4 shown in Fig. [51 When the sub- 
strate is strongly attractive (wet) with 7 = 0.4, the vol- 
ume fraction / increase linearly as the function of time 
step. In this case, the contact angle at the bottom 
substrate is determined from cos9 = 0.4/0.042 ~ 10, 
which gives complete wetting contact angle 9 = 0. The 
front velocity estimated from (|30|l gives Vg = 0.075 as 
D = 1/2 and e = 0.1, which nicely corresponds to the 
value Vg = 512/5500 ~ 0.093 estimated from the slope 




7000 



Time step r 

FIG. 6: The same as Fig. [3] but with a weakly wet substrate 
with 7i=-|-0.4, 0.1, 0.07 and 0.05 at the bottom for the low 
undercooling 7^=0.1 and at a medium temperature l/^o = 40. 
The standard deviations are indicated only for 7; = 0.07 and 
0.05 bu vetical lines. They are too small to be shown for 
7i = 0.1 and 0.4. 

1/5500 of the straight fine for 7; = 0.4 in Fig. M 

Therefore, in this regime, an Avrami plot ^ does not 
give any information about the microscopic picture of 
growth kinetics. It is also clear from Fig. [SI that when 
the substrate is very attractive (7/ = 0.4) the phase trans- 
formation is accomplished within a relatively short time 
compared to the weakly attractive case (7/ = 0.05). In 
the latter case, the bottom substrate acts to enhance the 
nucleation in the bulk, but the phase transformation still 
looks like homogeneous nucleation in Fig. [6l within the 
bulk. 



C. Homogenous-Heterogeneous coexistence regime 

When the undercooling is moderately high and the 
temperature is low, homogeneous nucleation and hetero- 
geneous nucleation can coexist if the surface field is at- 
tractive and strong. Figure [51 show a typical evolution 
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FIG. 7: The same as Fig[2]but with a strongly wet substrate with 7; — +0.4 at the bottom for the low undercooling rj — 0.1 
and at a medium temperature 1/Co = 40. Homogeneous nucleation is suppressed and nucleation occurs preferentially from the 
lower wet substrate. 



pattern when the surface field is 7; = 0.4. Since the tem- 
perature is low, the incubation time for homogeneous nu- 
cleation is long. Also, since the undercooling is medium, 
the growth velocity is slow. Then the heterogeneous nu- 
cleation and subsequent growth is not shadowed by the 
homogeneous nucleation. 

The corresponding evolution curve of volume fraction 
/ looks superficially similar to the KJMA prediction for 
the homogeneous nucleation and growth. However, the 
incubation times become apparently shorter as we in- 
crease the surface field 7; (Fig. [9]), which means that 
the phase transformation starts earlier because the in- 
cubation time for the heterogeneous nucleation becomes 
shorter from ([5]) and 

Figure (TU] shows the starting time to of phase transfor- 
mation as the function of the surface field 7;. Since we 
cannot use the KJMA formula ^ anymore to fit the sim- 
ulation data to deduce the incubation time Unc, we have 
used starting time to instead. The figure clearly indicates 
the expected decrease of starting time tg as we increase 
the surface field 7^. Therefore, the phase transformation 
can start earlier due to the presence of the wet substrate. 

The Avrami plots which correspond to Fig.[5]are shown 
in Fig. [TT] where we have used the starting time to esti- 
mated from Fig. [7] instead of the incubation time iinc in 
(m . Although the curves in Fig. [5] looks similar to the 
KJMA prediction, the Avrami plot in Fig. [TT] clearly dis- 
play a large deviation from the KJMA prediction of lin- 
ear relation in early stage of growth. When this surface 
field is strongly attractive (7; = 0.4), in particular, the 
Avrami plot consists of two linear segments with different 
Avrami exponents m. In the early stage of transforma- 
tion 771 ~ 1.06 while in later stage it becomes 777. ~ 2.94, 
which is very close to the ideal value 777 = 3 for two- 
dimensional isotropic growth predicted from ([7|). 

Therefore, nucleation at the early stage is dominated 
by heterogeneous continuous nucleation in dimensions 
or the site saturation in 1 dimension with Avrami ex- 
ponent 777 = 1 from ([8]). This contribution from hetero- 
geneous nucleation to the transformed volume / is soon 
suppressed due to the exhaustion of the nucleation site 
on the bottom substrate. At a later stage, heterogeneous 



nucleation is shadowed by homogeneous continuous nu- 
cleation which starts inside the bulk material since there 
is still much of the site to be use for homogeneous nucle- 
ation to start. Then the Avrami exponent in later stage 
becomes close to the ideal value 777 = 3 predicted from 
([7| for the d = 2 dimensional continuous nucleation. 

V. CONCLUSION 

In this study, we used the cell dynamics method fl^ 
to study the possible scenario of heterogeneous nucle- 
ation and growth of new phase in unified manner within 
the framework of the time-dependent Ginzburg-Landau 
(TDGL) model. We found that the formula for the vol- 
ume fraction of transformed volume due to Kolmogorov- 
Johnson-Mehl- Avrami (KJMA) theory and the so-called 
Avrami plot do not necessarily capture the essential fu- 
ture of the dynamics of heterogeneous nucleation. In par- 
ticular, when the heterogeneous nucleation is dominated, 
the Avrami plot is totally useless. Rather, a simple linear 
plot of the transformed volume / against time could give 
a linear relation between the transformed volume and 
time, and the front velocity directly gives the interfa- 
cial velocity of phase transformation. Therefore, a naive 
interpretation of Avrami plot when the supporting sub- 
strates are changed, needs careful examination [33l [34|. 

We have also predicted that, in general, heterogeneous 
nucleation and homogeneous nucleation could coexist. 
Then the simple Avrami plot would consist of two lin- 
ear segments with inflection. The Avrami exponents m 
of two segments would give the dominant mechanism of 
phase transformation. It is predicted that since the het- 
erogeneous nucleation could occur earlier, we might have 
a smaller Avrami exponent an at early stage and a larger 
Avrami exponent at a later stage. 

Since the TDGL model with the cell dynamics method 
is flexible and computationally simple, yet it has a direct 
connection to the thermodynamics of the system con- 
sidered, it can be used to study various scenarios of the 
dynamics of phase transformations. In fact it has already 
applied to those materials with a complex phase dia- 
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FIG. 8: The same as Fig. [2]but with a strongly wet substrate with 7/ — 0.4 at the bottom for a medium undercooling e = 0.2 



and at a lower temperature l/^o 



Initially the nucleation starts from the bottom substrate as the heterogeneous nucleation 



proceeds. Later, the homogenous nucleation and growth within the bulk starts, 
heterogeneous nucleation and growth coexist. 



Eventually, both the homogeneous and the 
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gram [35| . Although standard cellular automaton [36|, [STf 
and the lattice model could also be used to simu- 
late the complex dynamics of phase transformation, they 
have to introduce an artificial evolution processes algo- 
rithmically. On the other hand, the evolution in our cell 
dynamics method is driven by the free energy and the 
thermal noise, and it is completely free from artificial 
parameters. Therefore, our TDGL model with cell dy- 
namics method is more natural and easy to use in the 
study of various phase transformation scenarios, in par- 
ticular, of the time evolution including the incubation 
time without introducing many unknown and uncontrol- 
lable parameters. 



FIG. 9: The same as Fig. |3] but with a wet substrate with 
7i=-|-0.4, 0.1, 0.05 and 0.00 at the bottom for the medium un- 
dercooling 7^=0.2 and at a low temperature l/^o = 80. When 
the surface field yi increases, heterogeneous nucleation with 
a low energy barrier appears. Then nucleation starts earlier. 
The global shape of the transformation curve / does not dif- 
fer much from that of KJMA prediction for the homogeneous 
nucleation. However, the incubation time iinc becomes appar- 
ently shorter. Vertical lines indicate the standard deviations. 
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FIG. 10: The start time to ~ iinc of the phase transformation 
averaged over 50 samples estimated from the evolution of the 
volume fraction / as the function of the surface field 7; when 
e — 0.2 and 1/^q — 80. Vertical lines indicate the standard 
deviation. 
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FIG. 11: The Avrami plot corresponding to Fig. [S] Although 
the curves for the volume fraction / all look very similar to the 
curve for homogeneous nucleation in Fig[9l the Avrami plots 
clearly show the transition from one-dimensional growth at 
earlier time with a gentle slope with smaller exponent m ~ 
1.06 to the two-dimensional growth at later time with a steep 
slope with larger exponent m ~ 2.94 when the surface field is 
strongly attractive (71 = 0.4). 



